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ABSTRACT 



We have detected coherent osciUations, at multiple frequencies, in the line and 
^ ' continuum emission of the eclipsing dwarf nova V2051 Ophiuchi using the 10m Keck 

^\ I II telescope. Our own novel data acquisition system allowed us to obtain very fast 

OO . spectroscopy using a continuous readout of the CCD on the LRIS spectrograph. This 

'nJ" ' is the first time that dwarf nova oscillations are detected and resolved in the emission 

lines. The accretion disc is highly asymmetric with a stronger contribution from the 
blue-shifted side of the disc during our observations. The disc extends from close to 
the white dwarf out to the outer regions of the primary Roche lobe. 

Continuum oscillations at 56.12s and its first harmonic at 28.06 s are most likely 
• to originate on the surface of a spinning white dwarf with the fundamental period 

pinj corresponding to the spin period. Balmer and Helium emission lines oscillate with 

a period of 29.77s at a mean amplitude of 1.9%. The line kinematics as well as the 
^ • eclipse constraints indicate an origin in the accretion disc at a radius of 12 ± 2i?u,d. 

[ The amplitude of the emission line oscillation modulates (0-4%) at a period of 488s, 

. corresponding to the Kepler period at K—12Rt^d- This modulation is due to the beating 

between the white dwarf spin and the orbital motion in the disc. 

The observed emission line oscillations cannot be explained by a truncated disc as 
■ in the intermediate polars. The observations suggest a non-axisymmetric bulge in the 

^ I disc, orbiting at 12i?u,d, is required. The close correspondence between the location 

of the oscillations and the circularisation radius of the system suggests that stream 
overflow effects may be of relevance. 

Key words: accretion, accretion discs - novae, cataclysmic variables - stars: indi- 
vidual: V2051 Oph - stars: oscillations 
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1 INTRODUCTION 

Rapid variability is a tell-tale sign of accretion. The first 
coherent oscillation in a cataclysmic variable binary was re- 
ported by Walker as far back as 1956 (Walker 1956). The 
strictly periodic 71s brightness oscillations he found in DQ 
Her are now understood as being produced by a spinning 
magnetic white dwarf that is being fed by an accretion disc, 
and accretes the material along its magnetic field lines. Since 
the dynamical time scales of matter close to the white dwarf 
are of the order of tens of seconds, it was not until the devel- 
opment of high speed photon-counting photometers in the 
'70s that periodic signals were discovered in other systems. 



The dwarf nova oscillations (DNOs) are periodic brightness 
modulations discovered in dwarf nova systems during out- 
burst (Warner & Robinson 1972). The short time scales (2- 
40s) indicate an origin close to the white dwarf. In some 
systems, different periods were found at different times, and 
the period stability is considerably lower than the coher- 
ent oscillations in DQ Her. See Warner (1995) for review. 
Whereas the amplitudes in the optical are usually so low 
(< 1%) that the oscillations are only revealed after period 
analysis. Marsh & Horne (1998) detected coherent oscilla- 
tions at two periods in the HST light curves of OY Car 
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with an amplitude of up to 20% towards the end of a super- 
outburst. 

Reminiscent of the DNOs are the flood of QPOs re- 
ported in X-ray observations of accreting neutron star bi- 
naries (see van der Klis 2000 for review). The properties 
of these oscillations are in many ways very similar to those 
of DNOs, albeit on shorter time scales since these systems 
contain a much smaller compact object. It shows that mod- 
erately coherent oscillations appear to be a generic feature 
of accretion onto compact objects. 

V2051 Ophiuchi is a short period eclipsing dwarf nova 
system belonging to the SU UMa sub-class. Its status as 
a dwarf nova was contested (e.g. Warner & O'Donoghue 
1987) until two super-outbursts were observed in 1998 and 
1999 by amateur variable star networks. Spectroscopy re- 
veals double-peaked emission lines that exhibit a rotational 
disturbance effect during eclipse, corresponding to a pro- 
grade rotating accretion disc (Cook & Brunt 1983, Watts et 
al. 1986). Baptista et al. (1998) present a photometric model 
for the binary based on HST and ground based observations. 
Due to the unusual low-state the system was in during their 
observations, the orbital phases when bright spot and white 
dwarf are eclipsed by the mass donor star could be directly 
measured. 

Warner & Cropper (1983) concluded that most of the 
rapid variability (flickering) was produced in the inner disc 
regions rather than the bright spot or white dwarf. The only 
quasi periodic oscillation in this system was reported by 
Warner & O'Donoghue (1987), who detected a very short 
lived 42s oscillation during outburst. It was not classified as 
a DNO because of its low coherence, and the absence of a 
corresponding spike in the power spectrum at that period. 
We present high time resolution spectroscopy of the dwarf 
nova V2051 Oph shortly after a normal outburst, obtained 
using our own data collection hardware equipped on Keck 
II, and report on the discovery of emission line oscillations 
at the level of several percent. 

In the next Section we will describe the data acquisition 
system and novel reduction steps that were used. Section 3 
covers the analysis of the data and discusses the properties 
of the oscillations. We discuss likely models in Section 4 and 
summarise our conclusions in Section 5. 



2 OBSERVATIONS AND REDUCTION 

The observations were obtained as part of a 5 night cam- 
paign to study rapid variability of cataclysmic variables and 
X-ray binaries (O'Brien 2000). The optical data were taken 
using the Low Resolution Imaging Spectrograph (LRIS; Oke 
et al. 1995) on the 10-m Keck II telescope on Mauna Kea, 
Hawaii between UT 7:22 and 7:57 on July 4th 1998. The 
LRIS was used with a 5.2 arc-second slit masked with alu- 
minized mylar tape to form a square aperture. The 300/5000 
grating used has a mean dispersion of 2.55 A/pixel in the 
range 3600A- 9200A. We used a novel data acquisition sys- 
tem to obtain almost 28,000 2048 pixel spectra of V2051 
Oph, in the form of a continuous data stream. Each spec- 
trum has a mean integration time of 72.075 msecs and there 
was no dead-time between individual spectra. Each 2048 
pixel spectrum also contains a 25 pixel underscan region 
and a 75 pixel overscan region. The noise for a given pixel 



was calculated using a readout noise of 6.3 e~ and a gain of 
4.7 e~/ADU. Cosmic rays were rejected with a threshold of 
7-sigma from the de-biased frames. A master flat-field image 
was created by finding the median of 700 individual flat-field 
frames. This image showed no deviations above 0.3 % in all 
but 3 pixels. We decided that it was therefore not necessary 
to flat-field individual spectra. Calibration arc spectra were 
taken through the same slit aperture. 

Bacground spectra, including light from the sky but 
dominated by a constant electronic background, were taken 
at the beginning and end of target observing, so that long 
time scale variations could be detected. The mean and vari- 
able components of the background spectrum were found 
by creating a light curve for each pixel and extracting the 
mean and gradient of this light curve. These coefficients were 
then filtered in wavelength, with a running median filter of 
width 101 pixels. The gradient of the data was found to be 
almost zero. The thus derived background spectrum, which 
accounts for ^ 0.9% of the total flux, was subtracted from 
each target spectrum. 

The arc calibration was done by fltting a second order 
polynomial to 7 lines in a median spectra of exposures of Hg 
and Ar lamps. Arc spectra from the beginning and end of the 
exposures were used to take into account any drifts in the 
wavelength scale. The interpolated wavelength calibration 
was applied to each target spectrum using the MOLLyQ 
spectral analysis package. The spectrophotometric standard 
star, Feige 67 (Oke 1990) was observed using the identical 
setup. A low order polynomial fit was found to the median of 
all the individual flux star spectra and the individual target 
spectra were fiux calibrated using MOLLY. 

Individual time marks, accurate to 200ms, were placed 
after every other spectrum using the computer clock. In 
order to find the absolute times, secondary timestamps at 
known UT were placed at regular intervals throughout the 
observations. A calibrated UT radio clock was used for this 
purpose. With the help of these secondary timestamps, the 
times of each individual spectrum were calibrated, and we 
were thus able to achieve an absolute timing accuracy of 
close to 200ms. 

To monitor seeing conditions and telescope tracking, the 
auto-guider CCD output was stored. We used this output 
to measure the position and FWHM of the guide star image 
every ~2s during our observations. Gaussian fits to the stel- 
lar images reveal very stable tracking and seeing conditions 
throughout our observations. A conservative estimate of the 
seeing was derived from the mean FWHM of the guide star 
images. This suggests that seeing was better than 0.8 arc- 
seconds, equivalent to a seeing limited FWHM resolution of 
our spectra of 3.7 CCD pixels. 



3 DATA ANALYSIS 

Occasionally, not all pixels of a spectrum were successfully 
transferred to the acquisition PC. These incomplete spectra 
were rejected from the data stream. In order to increase the 

* MOLLY is a ID spectral analysis software package written and 
maintained by Dr. T. Marsh. Available for download at 
http : //www. astro . soton. ac .uk/ ^trm/ software .html 
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Figure 1. Top: trailed spectrogram of the whole data set after phase binning to 360 bins. The data are plotted on a logarithmic 
wavelength scale to provide a uniform velocity scale across the spectra. Gray scale denotes flux density in mjy. Bottom: A blow up of the 
eclipse of a few strong lines. Both Balmer and Helium lines show the classical pattern of the pro-gradely rotating accretion disc material 
being occulted by the companion star. Gaps occur at time of sky background measurements. 



signal to noise of the remaining spectra and make the data- 
set more manageable, individual spectra were binned to- 
gether in groups of 7 to give a final 3968 spectra with a time 
resolution of 0.5045 seconds. Since the dynamical timescale 
close to the white dwarf is of the order of a few seconds, 
this time resolution is still ample for our purposes. These 
spectra cover 40% of the binary orbit including an eclipse. 
The linear ephemeris of Baptista et al. (1998) was used to 
calculate orbital phases throughout this paper. Phase 0.0 
corresponds to inferior conjunction of the companion star 
i.e. mid-eclipse. 

Figure ^ shows a trailed spectrogram plot of the ex- 
tracted spectra obtained covering the orbital phases between 
0.69 and 1.08. The two gaps near phase 0.720 and 0.945 were 
introduced as the telescope was slewed to sky to allow robust 
sky subtraction as discussed above. Double-peaked emission 
lines are clearly visible, with the blue-shifted side of the lines 
eclipsing first, followed by the red-shifted peak, indicating 
the presence of a pro-grade rotating accretion disc. 

A brightening of the object to V—13.2 was reported 
by the VSNET^ network on July the 2nd (VSNET Alert 
#1934). The usual quiescent brightness level of this dwarf 
nova is around V=15.7. This suggests a normal outburst was 
in progress in early July. Unfortunately, no measurements 
were available for the preceding week that could confirm 
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the exact start of the outburst. At the time of our observa- 
tions, a V magnitude of 14.8 was reported, and 10 days after 
the observations the object was back to a quiescent level at 
V=15.7. It appears, therefore, we were observing V2051 Oph 
towards the later phase of decline from normal outburst. 



3.1 The average spectral properties 

Figure]^ (top) illustrates the average out of eclipse spectrum 
obtained by averaging all spectra between orbital phases 
0.69 and 0.92. Strong broad disc emission lines are super- 
posed on the flat continuum. The rapid drop in the spectral 
efficiency of the LRIS setup below 4000A made it difficult 
to achieve robust flux calibration for bluer wavelengths. The 
Balmer jump appears to be in emission, however care must 
be taken with the flux scale in this region. The Balmer lines 
have a full width half maximum of 2145 ± 20km s~^, and the 
Balmer decrement from Ha to HS is 1:0.65:0.50:0.41, sug- 
gesting a mix of optically thick and thin emission. The He- 
lium I lines are clearly double peaked and feature a slightly 
higher FWHM of 2245 ± 30 km s"^ The blue peak is con- 
sistently stronger during this phase interval. Curiously, only 
Hq exhibits a slightly stronger red peak. This blue asymme- 
try is also present in the average spectrum of Watts et al. 
(1986) that covers H/3 and H7, but not Ha. The Ha proflle 
in Warner & O'Donoghue (1987) again shows a stronger red 
peak. This curious reversal of the double peak asymmetry 
seems to be a persistent feature of V2051 Oph. The line 
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Figure 2. The average spectrum of V2051 Oph displays broad double-peaked Balmer and Hel emission lines. The top panel displays 
the out of eclipse spectrum obtained by averaging all spectra between orbital phases 0.69 and 0.92. The position of prominent lines as 
well as two telluric features red-wards of Ha are marked. During mid-eclipse (0.995 < <j> < 1.005), the continuum slope is slightly redder, 
and the emission lines only show the blue-shifted side of the disc as the red-shifted side is occulted from the observer (middle). The white 
dwarf contribution was isolated by subtracting the spectra on either side of white dwarf ingress and egress. It features a blue continuum 
and line emission from the accretion disc as the red-shifted part of the accretion disc is re-appearing (bottom). 



wings extend out to at least 3500 ± 200 km s~ for both the 
Balmer and Helium lines. 

Double Gaussians were fitted to the 6 strongest Balmer 
lines simultaneously. The double-peak separation was mea- 
sured to be «sep sin i = 1120 ± 60 km s~^, reflecting the 
radial velocities (±560 km s~^) of the gas on the outer edge 
of the accretion disc. If the velocity fleld of the disc is Kep- 
lerian and taking M^d = 0.78MQ,R^d = O.OlO4i?0 for the 
mass and radius of the white dwarf (Baptista et al. 1998), 
then the orbital radius Rkep corresponding to this velocity 
Vkep = Vsep /2 is givcu by : 

„ GM,,a f \fW00kms-^Y 
tlkcv = — n ~ 10 nr- cm 

This gives an outer disc radius of Rout = 3.3 x lO^^cm^ 
45±4J?u,tj. Although our estimated outer disc radius strictly 
exceeds the distance to LI {Rli = i2Ru,d), tidal interactions 
will make the disc velocities sub-Keplerian and keep the disc 
within the primary Roche lobe. Similarly, the emission line 
wings indicate that disc emission extends to ~ 3500 km s~^ 
corresponding to an inner disc radius of Ri„ ~ 8.4 x 10^cm= 
1.2 ± 0.2i?u,tj • In order to confirm that the emission lines are 
indeed formed in the disc, we measured the eclipse behaviour 
of the emission lines as a function of velocity. We extracted 
light curves for each pixel between ± 4250 km s~^ of II/3 and 
measured the eclipse half depth at ingress and egress (Fig- 
ure ^ . We see the blue-shifted side of the outer disc being 
eclipsed first followed by progressively smaller radii at larger 
velocities. The red-shifted disc emission is eclipsed later, 
with the outer edge disappearing just before mid-eclipse. 
The presence of a nearby Hel line at -1-3800 km s~^ distorts 



the measurements of the far red wing. Some uneclipsed light 
is present at low velocities. It is symmetric around zero ve- 
locity and has a FWHM of 1000 km s~^. Also, the average 
time of mid-eclipse occurs significantly earlier than white 
dwarf conjunction. This has also been reported in IP Pegasi 
during outburst (Steeghs 1999, Morales-Rueda et al. 2000). 
Analytical models for the velocity dependent eclipse phases 
of symmetric Keplerian discs are centered around phase 0.0 
(e.g. Young, Schneider & Shectman 1981). These models are 
clearly not able to reproduce the observed eclipse behaviour 
of V2051 Oph, which suggests a strong disc asymmetry is 
present. The analysis demonstrates that most of the line 
emission is indeed produced in a pro-gradely rotating disc 
that extends from close to the white dwarf to the outer re- 
gions of the primary Roche lobe. 

We isolated the mid-eclipse spectrum (Figure ^, mid- 
dle) by averaging all spectra between orbital phases 0.995 
and 1.005. As the relatively blue continuum from the in- 
ner disc is eclipsed, the cool companion star makes a larger 
fractional contribution, resulting in a slightly reddened con- 
tinuum slope. Weak emission lines with a very strong blue 
shifted emission component are present, indicating that the 
eclipse is not total and some (outer) disc emission is still 
visible around these phases. The fact that the blue-shifted 
side is much stronger than the red-shifted peak, points to 
an asymmetry in the outer disc that was already visible as a 
persistent asymmetry in the out-of-eclipse line profiles. Ei- 
ther there is excess line fiux produced in the blue-shifted 
side of the disc, or part of the red-shifted side is absorbed 
by a geometrically (and optically) thick structure. 
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velocity (km/s) 

Figure 3. An analysis of the eclipse of the H/3 line emission. 
Half depth ingress and egress was measured for each pixel in the 
line profile. Top panel plots the measured values (dots) as well as 
the time for mid-eclipse for each velocity defined as the middle 
between ingress and egress (line). For comparison, the white dwarf 
contact phases and conjunction are indicated by dashed lines. The 
dot-dashed line indicates the average phase of mid-eclipse for the 
emission lines across the profile. Bottom panel is the out of eclipse 
flux of each velocity as well as the minimum flux density. The 
instrumental resolution is indicated as a Gaussian profile with 
height 1. 
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Figure 4. Measuring the white dwarf contact phases. Top panel 
is a filtered version of the continuum light curve followed by its 
derivative. The fitted ingress and egress phases are marked by 
vertical lines. For comparison we also fitted the continuum slope 
between 4000 and 7000 A as a function of orbital phase (positive 
slope indicates a red spectrum). The occultation of the white 
dwarf is marked by a sudden steepening as its blue spectrum is 
occulted (second panel from bottom). The derivative of this slope 
is again used to measure the ingress and egress phases (dashed 
lines, bottom panels). 



3.2 The white dwarf 

In order to isolate the contribution from the white dwarf 
we subtracted the spectra from either side of white dwarf 
ingress and egress (Figure bottom). To correct, to first 
order, for the changing contribution from the disc, we per- 
formed a linear fit to each wavelength pixel on a short section 
just before and after ingress. This fit was used to extrapo- 
late the spectrum just before and after white dwarf ingress, 
and the white dwarf spectrum was obtained by taking the 
difference. The same procedure was applied at white dwarf 
egress, and the two derived white dwarf spectra, identical 
within the error bars, were averaged to derive our final white 
dwarf spectrum as plotted in Figure ^. Some emission line 
contribution from the disc is still present, although strong 
absorption cores of the Balmer lines are due to the white 
dwarf absorption line spectrum. 

The blue slope refiects the continuum contribution ex- 
pected from a hot white dwarf. Its relative contribution 
to the out-of-eclipse spectrum increases from ~ 15% at 
8OOOA up to ~ 30% at 4OOOA. Black-body fits to the con- 
tinuum shape indicates a white dwarf temperature around 
15000 ± 2000 K. Catalan et al. (1998) also derived a white 
dwarf temperature of 15000K from HST FOS data. 

In order to measure the contact phases of the white 
dwarf accurately, we applied the conventional method of us- 
ing the derivative of the continuum hght curve (e.g. Wood 



et al. 1986). The light curve was first filtered with a run- 
ning mean filter using a width of one third of the expected 
duration of the ingress and egress features. The numerical 
derivative was then calculated, and the ingress {(pwi) and 
egress (0me) are defined as the orbital phases of minimum 
and maximum derivative (Figure In order to determine 
these extrema, the derivative light curve was fitted with a 
Gaussian near the white dwarf features. The fitted centroids 
of these Gaussians are marked by vertical dashed lines. 

To exploit the available spectral information, we also 
fitted the continuum slope of each spectrum with a first or- 
der polynomial. The ingress and egress features are expected 
to show up prominently on the continuum slope as the blue 
white dwarf spectrum is eclipsed/reappears. Figure ^ plots 
the derived continuum slope after filtering using the same 
filter that was used for the continuum light curve. The slope 
is a fit to all continuum sections between 4000 and 68OOA, 
with all the emission line regions masked out. A positive 
slope corresponds to a spectrum that increases towards red 
wavelengths. Indeed, white dwarf ingress and egress are eas- 
ily identified as continuum slope changes. We performed a 
similar analysis to this light curve. The numerical derivative 
was calculated and the extrema were fitted with Gaussians. 

From the ingress and egress phase we derived the du- 
ration of the white dwarf eclipse; A<l)^d = 0me — 0mi and 
the time of mid-eclipse 00 = l/2((;/>u,e -I- <l)wi)- Our values 
(with a precision of ±2 x lO"**) as well as the values deter- 
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Table 1. WD contact phases 
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00 


lightcurve 


0.9635 


1.0315 


0.0680 


0.9975 


cont slope 


0.9642 


1.0313 


0.0671 


0.9978 


mean 


0.9639 


1.0314 


0.0676 


0.9976 


Baptista et al. 


0.9670 


1.0334 


0.0662 


1.00 



mined by Baptista et al. (1998) are listed in Table The two 
methods give compatible results, although the white dwarf 
ingress derived from the continuum slope is slightly later 
compared to the conventionally derived phase. We therefore 
averaged the results of the two methods to derive our final 
measurements. Our derived phases are systematically earlier 
compared to Baptista et al. (1998), indicating that the or- 
bital ephemeris used is off by ~12 seconds at the time of our 
observations. Our expected absolute time accuracy is better 
than 200ms, and this is therefore a significant difference. We 
did not attempt to construct a new orbital ephemeris but 
instead use our own measured white dwarf contact phases 
whenever needed as an internal calibration of the phasing. 
We also derive a larger value for A(/)„d, compared to Bap- 
tista et al. We remark that the system was observed just 
after outburst in our case, whereas Baptista et al. (1998) 
observed V2051 Oph during a rare low state. 

3.3 Time variability 

To study the rapid time variability of V2051 Oph we ex- 
tracted orbital light curves from the phase resolved spec- 
troscopy. All emission line regions were masked off to con- 
struct a continuum light curve spanning 4180-7960A(Figure 
^ top). We fitted a 3rd order spline to these continuum 
regions in order to determine the continuum shape and 
subtracted the fits from the data. From these continuum- 
subtracted spectra we then calculated emission line light 
curves. For the Balmer series we integrated all line fiux 
within ± 3500 km s~^ from HQ,H/3,H7,H(5,He and H8 using 
inverse variance weights to optimise the signal to noise in 
the resulting Balmer line fiux light curve (Figure middle). 
Similarly, we added the contributions of the six strongest 
Helium I lines marked in Figure ^ to construct a Helium 
light curve (Figure ^, bottom). 

Apart from the deep eclipse that is visible in all light 
curves, there is also variability on short time scales. No 
prominent pre-eclipse hump is visible, although the limited 
phase coverage prevents a more robust analysis of the pres- 
ence of bright spot emission in the light curves. The eclipses 
are not total, indicating some uneclipsed light is contributing 
to both the line and continuum emission (see also Figure |^. 
In particular, the emission line eclipses are highly asymmet- 
ric, and mid-eclipse occurs significantly before white dwarf 
conjunction. A strong disc asymmetry is apparently present, 
as expected earlier from the asymmetric emission line pro- 
files. A blow up of a small section of each light curve reveals a 
clear periodic signal in the Balmer line light curve with a pe- 
riod around 30s. As expected, the continuum is also affected 
by non-periodic fiickering, commonly observed in CVs. The 
measured white dwarf phases are marked as vertical dashed 
lines, and coincide with steep jumps in the continuum light 
curve. 



In order to search for (quasi)periodic oscillations, Scar- 
gle power spectra were computed from the lightcurves. The 
slow orbital variations due to the eclipse were subtracted us- 
ing a running mean filter with a length of 160 data points. 
The top panels in Figure ^ show the filtered light curve 
for the continuum and Balmer lines. Scargle power spec- 
tra between Hz and the Nyquist frequency at 1 Hz reveal 
clear oscillations at low frequencies. In the continuum, aperi- 
odic fiickering results in a large amount of power below 0.02 
Hz, i.e. on time scales of minutes - hours. The power rises 
steeply towards lower frequencies. A fit to the log(!/) ver- 
sus log(Power) plot of the unfiltered continuum light curve 
revealed a power-law index of -2.1 for this fiickering com- 
ponent, a common index observed in CVs (Bruch 1992). In 
Figure ^, this rising power at low frequencies is cut off by 
our filter. The peak at low frequencies in Figure ^ (top) is 
thus only a result of our filter and not an intrinsic peak in 
the power spectrum of the source. No significant power is 
present at the previously reported oscillation by Warner & 
O'Donoghue (1987) at 42.2s (0.024 Hz) in either the con- 
tinuum or the Balmer lines. However two other frequencies 
stand out near 0.035 Hz. In the Balmer line power spectrum, 
these same two periods are the sole significant frequencies 
in the power spectrum. The position of the two frequencies 
(/i=0. 03359 Hz, /2=0. 03564 Hz) as weU as their harmon- 
ics are marked by vertical lines in Figure ^. The continuum 
displays power at /i, /2 as well as /2/2 with similar power 
levels. Power at fi/2 is however absent. All of the peaks 
are well above the noise, as was confirmed by calculating an 
ensemble of power spectra while randomly jiggling the light 
curve points within their error bars. In order to check that 
the peak at /2/2 is not due to the increasing power at low 
frequencies, we used the fitted power law of the red noise 
(i.e. flickering) component mentioned above to measure the 
power in the peak relative to the red noise component (see 
also Figure ^) . The power at /2 /2 is at least a factor of ~3 
above the red noise component, so that its probability of 
occuring by chance is only ~ exp [—3] ~ 0.05. Furthermore, 
cutting the light curve in several segments and calculating 
power spectra for each, shows the presence of both peaks 
near /2 and /2 /2 in all segments. We are therefore confident 
that there are truly periodic components in the continuum 
light curve with the three frequencies /i, /2 and /2/2. 

The Balmer lines appear to oscillate solely at /i and 
/2. Interestingly, whereas /2 is stronger in the continuum, /i 
is strongest in the Balmer lines. We could already identify 
a clear oscillation in the raw Balmer light curve with fre- 
quency /i (l//i = 29.8 s. Figure ||). The amplitude of the 
oscillation varies considerably on a time scale of ~ 8 minutes, 
almost disappearing completely at some phases. This ampli- 
tude variation reminds us of the beating behaviour between 
two closely spaced periods. With /2 — /i = 0. 00205Hz or a 
period of 8.1 minutes, it seems the two frequencies /i and /2 
are indeed connected by the beat frequency fteat ~ f2 — fi 
corresponding to a beat period of Pbeat = 488s. 

The oscillations are coherent over our data segment. 
Since we cover only ~ 70 emission line oscillations, we can- 
not put strong constraints to the coherency of the periods 
(> 10^) in order to compare them with the typical coherence 
levels of DNOs (10* - 10®). 
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3.4 The origin of tlie oscillations 

The power spectra presented in the last section illustrated 
that there is a source driving the continuum and some line 
oscillations with a frequency of /2 . With equal power in the 
continuum at both /2 and its lower harmonic /2/2, it appears 
that this driving source has /2/2 as fundamental frequency, 
corresponding to a period of P2 = 2//2 = 56.12s. The emis- 
sion lines mostly respond at the frequency /i corresponding 
to a beat period of Pbeat ~ l/fteat = 488s between /2 and 
/i. Since all orbits in the binary are pro-grade, we expect 
only to see beats at lower frequencies {fi ~ f2 — fbeat) and 
not at f — f + fbeat if these periods are associated with or- 
bital motions. This suggests that Pbeat represents an orbital 
timescale associated with the emission line source. 

In order to extract some spectral information concern- 
ing the oscillation sources, we folded all our spectra on the 
two periods Pi = 2/fi = 59.53s and P2 = 2//2 = 56.12s. 
The phase zero point was set to the time of the first data 
point in our data. Before phase folding the data, we sub- 
tracted a running mean from all spectra to remove the non- 
oscillating background spectrum. This was achieved by ap- 
plying a running mean filter with a filter length of 120 spec- 
tra (60s) to the light curve at each wavelength pixel in the 
spectrum. Figure m plots the resulting trailed spectrograms 
for the two periods. On both periods, the continuum as well 
as the emission lines can be seen to oscillate. 

In the emission lines, the oscillation source is crossing 
from blue to red twice within each 59s, forming a tiger stripe 
pattern in the folded trailed spectrogram. Curiously, the line 
emission is mostly visible going from blue shift (~ -1000 
km s~^) to red-shift (~ 1000 km s~^). In Ha, the motions 
are strictly blue to red, the higher Balmer lines show a weak 
signature of a red to blue component as well. 



In order to measure the velocity amplitude of the tiger 
stripe pattern, we combined the dynamical data of several 
lines. For comparison. Figure |^ plots the Ha emission on its 
own, the summed contribution of the higher Balmer lines 
H/3-He as well as the sum of Ha-He. All Balmer lines share 
the same kinematics and phasing, though the Ha spectro- 
gram may be slightly different from the higher Balmer lines. 
In order to measure the velocity amplitude (K) of the spots, 
we measured the gradient {dv /dcp) of the tiger stripe as it is 
crossing zero velocity. If the spot traces a sinusoidal pattern 
as a function of phase, we then have K = tt- ^ at zero. This 

^ ' 2n d<p 

resulted in an estimated velocity amplitude of isT = 1210±80 
km s~^ from Ha, K = 1150 ± 60 km s~^ from the higher 
Balmer lines and K = 1207 ± 40 km s~^ from the sum of all 
the Balmer lines used. 

The azimuthal velocity of the orbit with Kepler period 
29.77s is 2800 km s~^, and clearly the emission line oscilla- 
tions can therefore not be attributed to disc material orbit- 
ing at this Kepler period. However, we already associated 
the line source with the beat period of 488s. If this period 
is a Keplerian orbit, it corresponds to a velocity of 1100 ± 
100 km s~^, very close to the observed velocity amplitude 
of the emitting material. 



3.5 Spectra of the oscillations 

We extracted oscillation spectra from the phase-resolved 
spectroscopy F{X,t). This was achieved by fitting sine func- 
tions with frequency / of the form; 

F{X, t) = A{X) sin (27r/t + ct>{X)) 

Here A{\) represents the oscillation spectrum and ^(A) 
the phase angle which was allowed to vary for each pixel 



© 2001 RAS, MNRAS 000, 



8 D.Steeghs et al. 





Ikyi 



0-04 0.06 
Frequency (Hz) 



lljL.IHi.ljiiLa..Lll.l.,.,..inUUll ,l...l.^.ui.L.l.U.J...I.IUkjli..........i.......,Ull 



^ikLiVj-AJ AJJ-L JiiLjjLL..li...J .UU.. .JlllA.MLIllL.lt 



0.4 0.6 

Frequency (Hz) 




fViiWA..i«.J.-».t.-i..t 



0.02 0.04 0.06 0.08 

Frequency (Hz) 

i.li.,.]J^.....j]i;.i...Ji,l...i,ij.l.. ... i...i.alluiLI.Ul.ll i. - .■.■.....i...J.J.J ..,[,.i.i.iUll.i..iLi,U. lL..ll.i..,.U.IL..,. 



0.2 0.4 0.6 0.8 

Frequency (Hz) 

Figure 6. Power spectra calculated from the continuum (top) and Balmer light curves. The top panel in each plot is the light curve 
after filtering with a running mean. For plotting purposes, the data are binned by 4. Main panels are Scargle power spectra calculated 
from the un-binned filtered data. A zoom in of the region between 0.01 and 0.08 Hz illustrates the detected oscillations. Dashed lines 
denote the position of the two frequencies and their harmonics are marked by dot-dashed lines. For the continuum, a power-law fit to 
the red noise component is indicated by the dashed curve. 



since. Figure |9| plots the fitted oscillation spectra obtained 
by fitting simultaneously to the frequencies /i and /2. We 
only used the data between phases 0.7 and 0.85 in order 
to avoid any problems with phase changes during eclipse. 
Since the signal to noise in the lines is much better com- 
pared to the continuum, we rebinned the spectra to obtain 
a minimum signal-to-noise of 50 in each bin. This allows 
us to keep narrow bins across the lines and wider ones for 
the continuum. Both continuum and lines oscillate at these 
periods as expected from the power spectra, with the lines 
achieving amplitudes of up to 5 times those of the contin- 
uum. Line emission modulated at /i appears to come from 
slightly higher velocities compared to the emission from /2, 
producing a clear double peaked emission line profile in the 
Balmer series and Helium lines. Since the phase angle is a 
free parameter for each pixel the amplitudes contain a pos- 



itive bias. However, the continuum phases are constant as 
a function of wavelength, and continuum amplitudes can 
thus be reliably obtained from the orbital continuum light 
curves. In the next Section, we will look at the variation of 
the oscillation amplitudes as a function of orbital phase. 

3.6 Eclipse constraints on the location of the 
oscillations 

The eclipse provides valuable constraints on the possible lo- 
cations of the oscillations. From the orbital light curves it 
is already clear that the line oscillations persist well into 
the eclipse. We constructed two Balmer emission line light 
curves to separate the blue-shifted and red-shifted emission. 
From the continuum-subtracted spectra we integrated all 
line flux between -3500 and km s~^ of Hoi — H%. Similarly, 
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line kinematics across the spot phase. 



we constructed a light curve for all the red-shifted line emis- 
sion between and -1-3500 km s~^. These hght curves were 
then filtered again by a running mean filter of 80 points. 
We then fitted sinusoids to sections of the hght curves to 
determine the amplitudes and phase of the oscillation as 
a function of orbital phase. Since the two frequencies are 
so close to each other, it was not possible to reliably fit the 
amplitudes of two frequencies simultaneously. Instead we fit- 
ted the amplitude of the strongest line oscillation at Pi/2 
(29.8s). Figure shows the derived amplitudes for the blue 
and red hght curves as well as a re-binned version of the 
two light curves to illustrate the modulations in amplitude 
directly. The amplitudes peak at a level of ~ 4% and mod- 



ulate regularly with the beat period of ~ 8 minutes. There 
also appears to be a general trend of a gradual drop in the 
oscillation amplitudes across the data, indicating the oscilla- 
tion has a limited lifetime. The oscillation in the blue wings 
drops down around orbital phase 0.93, and only slowly re- 
covers after mid-eclipse. On the other hand the oscillations 
persist in the red wing until mid-eclipse, well past the white 
dwarf ingress at phase 0.963. This indicates that the site 
of the emission line oscillations cannot be confined to the 
white dwarf, but needs to be located at a considerable dis- 
tance from it. Only red-shifted emission from the oscillation 
is visible during the early parts of the eclipse. This is in ac- 
cordance with an origin as a pro-grade rotation in the disc, 
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Figure 10. The Balmer oscillation amplitudes separated into blue-shifted (left) and red-shifted emission. Top panels show the unfiltered 
light curves, followed by a filtered and binned version to highlight the oscillations. The amplitude variation can be directly identified in 
the two light curves. Sliding sine fits to the data were used to obtain a light curve of the oscillation amplitude. The fitted amplitudes 
and phases constitute the next two panels. Final panel plots the RMS of the fit residuals after the sine-fit is subtracted from the data. 
Dashed lines mark white dwarf ingress and ogress. 



since the blue-shifted part is already eclipsed, but the red- 
shifted part is still visible. This also fits in well with the ob- 
served kinematics of the emission lines (Section 3.4, Figures 
^ ^and|^. The corresponding Kepler orbit is at a distance 
of 12±2i?„d. The phases of the oscillation modulate in sync 
with the amplitude changes across each beat cycle with a 
variation of ~ 100° . The blue side consistently advances the 
red by ~ 90° . If the emission line is indeed produced by ma- 
terial in the disc, phase changes are expected as the emitting 
material orbits around the white dwarf During eclipse both 
sides show a large change in the phasing of the oscillation 
as part of the oscillating source is eclipsed, with the two 
sides being out of phase by roughly ~ 280° at the end of our 
data. We also calculated the RMS of the fit residuals after 
the sine fit was subtracted from the data. The RMS light 
curve shows a short interval of increased variability, near 
white dwarf ingress for the blue side, and near egress for the 
red-shifted side. It appears to originate from the small area 
of the inner disc that is still visible, in between the Roche 
lobe of the companion and the white dwarf. The variability 
is not periodic and we do not seem to see similar variability 
from the side of the disc on the other side from the white 
dwarf (i.e. at egress for blue and ingress for red). 

The amplitude/phase measurements for the continuum 
oscillation were more complicated to extract. First of all, as 
can be seen from the oscillation spectra, the amplitude of the 
continuum oscillation is typically about a factor 5 weaker, 
and there is also a lot of aperiodic variability in the con- 
tinuum superposed on the oscillation. Since the phasing of 
the oscillation appears to be wavelength independent (Fig- 
ure bl) , we do not need to split the continuum emission up 



Table 2. Fitted oscillation amplitudes 

amplitude range mean amp 

Balmer lines 0.01-0.26 mjy (0.0-3.7%) 0.11 mjy (1.9%) 

continuum 0.00-0.04 mjy (0.0-2.%) 0.022 mjy (0.8%) 



into wavelength bands, and can use the orbital continuum 
light curve we derived before. The continuum light curve 
was first filtered using a running mean to remove the slower 
aperiodic variations to first order. Sine functions were then 
fitted at the continuum oscillation period of P2 = 56.12s. 
Unfortunately, the continuum was clearly affected by other 
aperiodic variability as the fits were poor with high values 
(~ 3-5) and large RMS residuals. The oscillation amplitude 
and phase varies considerably both during eclipse as well as 
out of eclipse. During eclipse the amplitudes drop and are 
compatible with zero. However, no sharp eclipse of the os- 
cillation was apparent that could firmly establish an origin 
on the white dwarf itself The amplitude light curve of the 
continuum oscillation is therefore not accurate enough to 
distinguish between a location on the white dwarf itself or 
close to it in the disc. 

Table ^ lists the oscillation amplitude ranges obtained 
from these sine fits to the Balmer lines and the continuum. 
Although the continuum oscillation is much weaker com- 
pared to the lines, both reach maximum relative amplitudes 
of a few percent. We have obtained some strong constraints 
to the origin and location of the continuum and emission 
line oscillations. In the next section we review several possi- 
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Figure 9. The spectrum of the oscillation obtained by fitting 
sine-functions to the data for each pixel. Top two panels are spec- 
trum of the oscillation with frequency /i and the phase angle, the 
bottom two panels correspond to /2. The oscillation spectra are 
rebinned to achieve a minimum signal to noise in each spectral 
bin of 50 in order to exploit the better s/n in the lines compared 
to the continuum. 

ble models for the observed oscillations, and see if they are 
able to satisfy the above constraints. 

4 INTERPRETATION 

Before reviewing possible explanations for the observed os- 
cillations, let us quickly summarise the main features that 
a model must satisfy in order to explain the observations. 
In the continuum we observe oscillations at 56.12s and its 
harmonic at 28.06s, that must originate close to the white 
dwarf or on the compact object itself. The emission lines os- 
cillate mostly at 29.77s, and require a source at a distance of 
~12f?™tj, corresponding to an orbital period of 488s, which 
is also the beat period between 29.77s and 28.06s. The emis- 
sion line oscillation amplitudes modulates regularly on the 
beat period of 488s and emission is only visible going from 
maximum blue-shift to red-shift, but not vice- versa. 

4.1 Pure disc oscillations? 

Over the years various models have been proposed that at- 
tribute the DNOs to oscillations or modes in the disc it- 
self. Possible mechanisms are non-radial pulsations (e.g. van 
Horn, Wesemael & Winget 1980), axisymmetric radial pul- 
sations (Kato 1978), vertical oscillations or resonant oscil- 
lations due to shock waves (Molteni et al. 1996). In most 
cases, the oscillation period is to first order the local Kep- 
lerian period. We can then obviously connect each observed 
frequency to a particular radius in the disc which gives orbits 



between 1 and 3 white dwarf radii out for the observed os- 
cillations. For example, P2 corresponds to a Kepler orbit at 
R = 2.8ii,i,d, and its harmonic is at i? = LSR^nd- However, 
there is no particular reason why the disc would oscillate at 
these distinct radii, and not at some other, and no explana- 
tion of how it manages to do so coherently. Also, the pres- 
ence of the harmonic of /2 but not /i is not explained, nor 
is the connection between /i, /2 and the observed beating 
behaviour with a period of 488s. More seriously, we showed 
in Section 3 that the emission line oscillation originates at 
12f?™tj, but oscillates with a period of 29.77s, which is not 
the local Kepler period. Although we have established that 
the emission line oscillation requires an origin in the disc, 
which is an important result in itself, a model relying on 
disc modes only seems unsatisfactory. 

4.2 White dwarf pulsations 

The emission line oscillation cannot originate from the white 
dwarf itself, but the continuum oscillations at P2 and P2/2 
are indeed expected to be related to the white dwarf. One 
possible mechanism that can produce periodicities apart 
from a fixed structure on a rotating white dwarf, is oscil- 
lations of the white dwarf. 

Papaloizou & Pringle (1977) discussed non-radial oscil- 
lations of rotating stars and concluded that oscillations in 
rotating white dwarfs can indeed produce the oscillations at 
periods corresponding to DNOs. The periods of these modes 
with order m is ~ Papin/m, with Papin the spin period of the 
white dwarf, and in principle any number of orders can be 
superposed at any given time. It is thus not obvious how one 
mode will stand out compared to the others and produce a 
single, coherent spike in the power spectrum. 

The detected power at P2 and P2/2 could thus, in prin- 
ciple, be interpreted as the incarnation of two orders of 
a non-radial oscillation on the white dwarf. However, two 
equatorial spots on a rotating white dwarf with spin period 
P2 could explain the observations just as well. The latter 
would also explain the observed beating behaviour between 
the spinning white dwarf and orbiting disc material more 
naturally. One way to distinguish between these two models 
would be the detection of other periods, but not P2, in this 
system. In the case of white dwarf oscillations it is easy to 
produce power at different orders and modes, however the 
white dwarf spin period cannot be changed on short time 
scales and should leave a persistent signature. 

Another class of models that relies on the origin of the 
DNOs close to the white dwarf are differentially rotating 
surface layers above the white dwarf (Warner 1995), or in- 
teraction between the disc and the boundary layer around 
the white dwarf (e.g. Popham 1999). The main concern with 
these models, when trying to understand the behaviour of 
V2051 Oph, is that they can in principle account for the os- 
cillations that originate close to the white dwarf (/2,/2/2), 
but not the emission lines oscillation (/i). 

4.3 An intermediate polar scenario 

Emission line oscillations have been seen so far only in the 
CV sub-class of intermediate polars such as DQ Her (Martell 
et al. 1995) and RXJ0558 (Harlaftis et al. 1999). In these 
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systems, the magnetic white dwarf accretes via two hot spots 
at its magnetic poles as it spins rapidly. The disc is truncated 
at the radius where magnetic pressure is comparable to the 
gas ram pressure, and the gas flow makes a transition from 
Keplerian motion in the plane to a motion along field lines. 
The observed line oscillations are produced by disc material 
that is illuminated by the hot spots on the white dwarf as 
they sweep by. 

However, if we try to interpret the observed oscillations 
in V2051 Oph in terms of this model, it fails for several rea- 
sons. We have established that the emission lines are formed 
in the disc which extends from very close to the white dwarf 
out to the outer edges of the primary Roche lobe. The emis- 
sion line oscillation is produced in the disc at ~ 12Rwd, 
but the disc is clearly not truncated there. In this model, 
P2 is the spin period of the white dwarf, which leads to a 
co-rotation radius of Rco = 2.8Rwd- Apart from the fact 
that the disc apparently extends down to the white dwarf, 
the system would be a propeller and any material thread- 
ing at _R ~ 12Rwd would be propelled out. Finally, in the 
intermediate polars we observe the emission lines as they 
move from red-shift to blue shift and not vice versa as in 
the case of V2051 Oph. If we would interpret the observed 
line oscillation kinematics (Figures in terms of material 
moving along magnetic field lines rather than in the disc, we 
wouldn't see the blue wing oscillations being eclipsed first 
followed by the red wing at mid-eclipse (Figure hol). The 
intermediate polar scenario therefore fails in the sense that 
the emission line oscillations cannot originate from the edge 
of a a truncated disc as in DQ Her and RXJ0558 or from 
a magnetic curtain. However, the system may well contain 
a weakly magnetised white dwarf with two active magnetic 
poles that is responsible for the oscillations at /2 and /2/2. 
However, this cannot account for the emission line oscilla- 
tions at 29.77s. 

4.4 A spinning white dwarf illuminating the disc 

We have reviewed some possible models that could explain 
the observed set of oscillations in V2051 Oph. The contin- 
uum oscillations at /2 and /2 /2 are most naturally explained 
by a spinning white dwarf with Pspin = P2 containing two 
hot spots, that produce a modulation at the spin period and 
its harmonic. The origin of these hot-spots may be due to a 
weak magnetic feld of the primary, which is low enough such 
that the disc can still extend to close to the white dwarf. The 
radiation from these two spots illuminates the disc, and is 
reprocessed as continuum and line emission at a radius of 
~ 12 ± 2R^d- The corresponding emission oscillates at the 
beat frequency between the spin of the white dwarf and the 
local Kepler period of 488s producing line oscillations at 
29.8s. 

We noted before that we only see line emission going 
from blue shift to red shift and not vice-versa. This may 
be caused by the presence of material close to the white 
dwarf above the disc plane, that absorbs the line emission 
when it is on the far side of the disc from the point of view 
of the observer, but not when it is in front. One candidate 
for this is the uneclipsed light. Alternatively, the emission 
of orbiting material could be beamed an-isotropically, only 
producing line emission away from the white dwarf. Orbiting 
material would then only be visible going along the front 
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Figure 11. A top view of V2051 Oph with the mass donor star 
on the right. The ballistic stream trajectory is plotted down to its 
point of closest approach to the white dwarf. The dot-dashed line 
is the site of the emission line oscillations and corresponds to the 
circularisation radius. The arrows on the outer disc denote the lo- 
cal direction of the Kepler flow. The arc on the outside illustrates 
the orbital phases during which the system was observed (time 
increasing clockwise). 

side. Perhaps the inner side is blasted by the white dwarf 
and is too hot to produce optical line emission, whereas the 
outer side is at low enough temperatures to emit Balmer 
and Helium emission. 

What remains to be explained is why does this process- 
ing occur at this specific radius, and why does it modulate 
its amplitude regularly? Since the orbiting disc material at 
this radius is beating with the white dwarf spin, it needs 
to be a localised region, both in terms of radius as well as 
azimuth, that is mostly responding for it to produce a coher- 
ent oscillation. A localised 'bulge' in the disc at the required 
radius would orbit the disc every 488s and its oscillation 
amplitude would modulate because of the intervening gas 
that prevents us from viewing the bulge when its on the far 
side. This obscuring material that was required to explain 
the blue-to-red pattern in the emission lines, thus also ex- 
plains the regular modulation of the amplitude every 488s. 
Popham (1999) also invoked a non-axisymmetric bulge in 
the disc to produce DNOs. In his model, the bulge is formed 
at the transition radius between the disc and the boundary 
layer surrounding the white dwarf. In our case, the bulge 
is clearly too far from the white dwarf to be due to such a 
boundary layer-disc interaction. We are now in the familiar, 
and unfortunate, position again that one needs to invoke 
some unknown mechanism in order to produce a structure 
in the disc at this specific radius. 

Figure ^ plots a schematic top view of the binary, with 
the site of the line oscillation marked as a circle. We calcu- 
lated a few characteristic radii for V2051 Oph, using the bi- 
nary parameters of Baptista et al. (1998). For example, we 
couldn't connect this radius with a tidal resonance. How- 
ever, the circularisation radius for this system is (Warner 
1995); 

Rc^rc = 0.0859g~° '*^^a = 7.8 x lO^m ~ 11 ± 2R^d 

This close correspondence between the emission line source 
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and the circularisation radius suggests that gas stream ef- 
fects may be of relevance. Lubow (1989) explored the possi- 
bility of stream overflow, and in the case of no disc interac- 
tion streamlines typically re-impact in the orbital plane at 
the radius of closest approach; 

-Rmin = 0.0488g~'' '"'''a = 4.9 x let'em ~ 7 ± IR^a 

The exact nature of the stream-disc interaction is not well 
understood, but we may expect some stream-disc effects be- 
tween these two radii. The attractive feature of such a stream 
overflow would be the natural way in which a vertically ex- 
tended structure would be present at roughly the right ra- 
dius to produce the reprocessed emission. We also noted in 
section 3.1 that near mid-eclipse the blue shifted contribu- 
tion from the disc is persistently stronger compared to the 
red-shifted emission. Stream overflow would block part of 
the red-shifted disc emission from view, and this asymme- 
try could thus be easily accommodated. The reversal of this 
asymmetry in the Ha remains a puzzle, however. 

It is clearly not trivial to produce a localised bulge in 
the disc at this stream impact location, that is able to orbit 
coherently for quite some time while stream overflow is oc- 
curing. Observations of this kind spanning longer intervals, 
could shed some light on the coherency of such structures. 
A more thorough simulation of such a scenario is also re- 
quired in order to test if such a geometry is a viable way 
of producing the observed oscillations (and perhaps other 
DNOs). This is beyond the scope of this paper, but such 
a scenario does appear to be able to accommodate all the 
features of the observed oscillations. Since we observe the 
system towards the end of an outburst, the radial density 
distribution is not expected to be that of a steady-state disc. 
Disc instability calculations usually show a broad peak in the 
density at roughly the circularisation radius towards the end 
of the outburst. Such an 'empty' disc would be in particular 
susceptible to stream overflow, and again the obvious region 
for stream-disc interaction would be near the circularisation 
radius where the density peaks and the stream is flowing 
towards the orbital plane again. 



5 CONCLUSIONS 

We detected coherent continuum and emission line oscilla- 
tions in the dwarf nova V2051 Oph on decline from a normal 
outburst. Accretion disc emission extends from very close to 
the white dwarf out to the outer parts of the primary Roche 
lobe. The disc emission lines display a persistent blue to 
red asymmetry, with the blue peak being stronger in all the 
Helium lines as well as the Balmer lines, except Ha which 
has a stronger red peak. The eclipse light curves are also 
highly asymmetric which suggests that the blue side of the 
disc makes a larger contribution to the emission compared 
to the red-shifted side. 

The continuum oscillations are most likely to originate 
on the surface of a spinning white dwarf with spin period 
56.12s and temperature around 15000K. The amplitude of 
the oscillation in the continuum varies between and 4%, 
and disappears during white dwarf eclipse. The Balmer and 
Helium I emission lines oscillate strongly at a period of 29.8s. 
The line kinematics as well as the eclipse constraints indicate 
these to come from a non-axisymmetric bulge in the disc at 



a radius of 12 ± 2Rwd- The corresponding Kepler orbit has 
a period of 488s, and corresponds well with the observed 
amplitude variations (0-4%) on this period. The oscillating 
line emission is observed to go from maximum blue-shift to 
maximum red-shift every 29.8s, but not vice versa. This, to- 
gether with the regular modulation of the oscillation ampli- 
tudes can be explained by intervening gas above the orbital 
plane, close to the white dwarf. This is supported by an un- 
eclipsed component in both the continuum and lines. In the 
Balmer lines this component is centered on zero velocity and 
has a FWHM of 1000 km s~^. Alternatively, the emission is 
beamed away from the white dwarf. 

The close correspondence between the location of the 
oscillations and the circularisation radius of the system as 
well as the disc asymmetry may indicate the relevance of 
stream overflow to the presence of vertically extended bulges 
in the disc. More observations of this kind as well as more 
detailed simulations of such a scenario are required to con- 
firm its feasibility in the light of producing DNOs. Predic- 
tions of our interpretation are the persistent presence of 
56.12s/28.06s oscillations from the white dwarf. This period 
should be detected at other epochs, most likely in the UV. 
The UV also provides opportunities to measure the spin of 
the white dwarf directly using the rotational broadening of 
white dwarf absorption lines (Sion 1999). With a spin pe- 
riod of 56.12s, the expected surface velocity of the spinning 
white dwarf works out to be: 

v^ud = f2 X = 'ii^/Pi X R^ud ~ 810 km 

Secondly, we expect the period of the emission line oscil- 
lations to change across the outburst cycle as the location 
and extent of stream overflow and re-impact will depend on 
the radial density distribution of the disc. For example, if 
we wish to interpret the reported 42.2s oscillation (Warner 
& O'Donoghue 1987) in this scenario, it can be identifled as 
the beat between the white dwarf spin and a Keplerian orbit 
at ~ 7Rwd- 

The data presented here demonstrate the potential of 
high speed spectroscopy of CVs using large aperture tele- 
scopes. We detected emission line oscillations in a dwarf 
nova for the first time, and located their origin in the ac- 
cretion disc. A better understanding of the elusive DNOs 
would benefit from a high time resolution study of (eclips- 
ing) dwarf novae across their outburst cycle. 
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